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Zj ' Abstract 

We consider Friedmann-Lemaitre-Robertson- Walker flat cosmological models in 
bJT)| the framework of general Jordan frame scalar-tensor theories of gravity with arbi- 

trary coupling function and potential. For the era when the cosmological energy 
density of the scalar potential dominates over the energy density of ordinary mat- 
ter, we use a nonlinear approximation of the decoupled scalar field equation for the 
— |- regime close to the so-called limit of general relativity where the local weak field 

CN ■ constraints are satisfied. We give the solutions in cosmological time with a particu- 

lar attention to the classes of models asymptotically approaching general relativity. 
The latter can be subsumed under two types: (i) exponential convergence, and 
(ii) damped oscillations around general relativity. As an illustration we present an 
example of oscillating dark energy. 

^ 1 Introduction 

In our recent paper pQ we considered the Friedmann-Lemaitre-Robertson- Walker (FLRW) 
flat cosmological models in the framework of general scalar-tensor theories of gravity 
(STG) with arbitrary coupling functions, set in the Jordan frame, in the era when the en- 
ergy density of the scalar potential dominates over the energy density of ordinary matter. 
We used a nonlinear approximation of the decoupled scalar field equation and analyzed 
the corresponding phase portraits in the neighbourhood of the point in the phase space 
which corresponds to general relativity (GR). We argued that the topology of phase por- 
traits of approximate nonlinear scalar field equation is similar to the topology of phase 
portraits of the exact equation. In the present paper we extend the study to the evolution 
in the cosmological time. 
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According to contemporary cosmological precision data, the expansion of the universe 
has been accelerating for the last few billion years. The phenomenon can be accom- 
modated in a specific FLRW model with a suitable cosmological constant A which is 
equivalent to an additional matter component with an equation of state (EoS) p = wp, 
w = — 1. However, several recent data analyses indicate that the dark energy (effective) 
barotropic index w is actually changing in the cosmological time [21 [3]. In particular, 
there may have been periods with dark energy in the phantom regime, w < — 1, which 
has inspired investigations of cosmological models allowing to cross the phantom divide 
value w = — 1 jl]. Alternatively, according to some authors [5j [6] the acceleration may 
already have passed the peak and is currently slowing down, while others have looked for 
the oscillating patterns in the dark energy which does not cross the phantom divide [TJ. 
Related theoretical models have been also proposed [8]. 

In this paper we present approximate solutions of the STG FLRW equations in the 
era when the energy density of the scalar potential dominates over the energy density of 
ordinary matter and demonstrate that depending on the choice of the model parameters, 
they include crossing the phantom divide as well as oscillating dark energy. The plan of 
the paper is the following. We begin in Sec. 2 with the general equations of STG FLRW 
cosmological models and their nonlinear approximation introduced in our previous paper 
[1] . Sec. 3 presents the solutions of approximate decoupled scalar field equation and gives 
an additional hint that the same types of solutions occur also in the full theory. An 
example of oscillating effective barotropic index of the dark energy is given in Sec. 4. 
Finally, Sec. 5 provides concluding remarks. 



2 Equations 



We consider a general scalar-tensor gravity in the Jordan frame governed by the action 
functional 



S 



1 

2^2 



d Xy 



-g 



VR{g) - ^-W P W P ^ - 2kV(¥) 



(1) 



Here oj(*f>) is a coupling function (we assume 2u(t) + 3 > to avoid ghosts in the Einstein 
frame, see e.g. Ref. [9]) and V(\&) > is a potential, V M denotes the covariant derivative 
with respect to the metric g^ u , S m is the matter action, and k 2 is the non-variable part 
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of the effective gravitational constant In order to keep the latter positive we assume 
that < * < oo. 

The field equations for the flat FLRW line element 



ds 2 = -dt 2 + a(tf (dr 2 + r 2 {dd 2 + sin 2 6dip 2 )) 
and barotropic fluid (p = wp) read 
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where H = a/a. 

We can express the Hubble parameter if as a function of ^ by solving the Friedmann 
equation ([3]) algebraically, 



J=4J W f) + »)-^ + ^. (6) 
2^ V ; 12^ 2 3# v ; 

Upon introducing the notation 

and substituting Eq. (jHJ) in Eq. (jHJ), we get a decoupled equation for the scalar field 

* = (^ + ^(*»( 2w (*) + 3 ))* 2 + ^MT3( 1 - 3w )" 

± ^ + + i^wffl * + 3 ■ (8) 

In the limit ( 2 ^(^)+3) ~ ^ 0, ^ 7^ the system faces a spacetime curvature singularity, 

since H diverges, and likewise behaves At first the limit (a) ( 2 u(4)+3) ~ >• 0, (b) \l/ — >■ 
seems only slightly less mathematically precarious for the equations are left just indeter- 
minate (contain terms jj) . Yet the latter situation is of particular physical importance, 
as the experiments in the Solar System (where matter density dominates over the scalar 
potential), i.e. the limits of observed values of the parametrized post-Newtonian (PPN) 
parameters and the time variation of the gravitational constant |10j . 

k 2 2w + 4 
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2 dw 
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u + 2 
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suggest the present cosmological background value of the scalar field to be very close to 
the limit (a)-(b). As in this limit the STG PPN parameters coincide with those of general 
relativity we may tentatively call (a)-(b) 'the limit of general relativity' 

Let us define ^ by 2cJ (^ ) + 3 = 0- l n our previous papers [HUH [12] we studied the limit 
(a)-(b) with the simplifying assumptions (c) A* = A(^/^) ^ and (d) is differentiable 
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at which enabled to Taylor expand the functions in Eq. (|SJ and find analytic solutions 
in the phase space for the resulting approximate equation. The outcome was that the 
solutions are well behaved in this limit, motivating the inclusion of = \I/ = 0) as 
a boundary point to the open domain of definition of Eq. (jSJ. Moreover, it was possible 
to identify a wide class of STGs where the FRLW cosmological dynamics spontaneously 
draws the scalar field to this limit, i.e. into agreement with current local weak field 
observations in the Solar System. 

The aim of the present paper is to complement our previous study [1J of the cosmologi- 
cal epoch when the energy density of the scalar potential dominates over the energy density 
of the ordinary matter by considering the limit (a)-(d) and finding analytic solutions not 
in phase space variables but explicitly in time, in order to facilitate the assessment of 
their cosmological viability. Let us begin by introducing ^(t) = ^ + x(t), and assuming 
that x(t), x(t) are both first order small. Then 

i i 

- + A*x(t) + ... « A*x(t) . (13) 



2w(*(t)) + 3 2u;(#, 
By denoting W* = W(^*) and 



C 1 = ± x l 3 ^^, C 2 = AW*, (14) 
the decoupled scalar field equation (JSJ) with p = is approximated by 

x + Ci x - C 2 x = — . (15) 
2x 

We will investigate its solutions in the next section. 
Analogously we may express 

rr Cl 1 ■ , 1 (CI C 2 \ 1 X 2 , . 
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A necessary condition for crossing the phantom divide w e g = — 1 is vanishing of the second 
term in Eq. (ITS]) . This occurs if A equals to a solution of the corresponding quadratic 
equation 
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Here we can see a condition on the constants A*, Ci, C*2 which characterize the theory: 
must be real numbers, i.e. D > 0. Note that on the plane (x, x) the phantom divide 
(TT9]) consists of two straight lines crossing at the origin x = 0, x = 0. 
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The PPN parameters ffTUj) . (TTTT) and ^ f lT2|) are approximated by 
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Strictly speaking, for constraining cosmological models using PPN type observations re- 
quires solutions with p ^ 0. In the following we will consider cosmological solutions with 
p = 0, V(^) 7^ for the late time universe, but since the matter density in the Solar Sys- 
tem exceeds considerably the cosmological matter density local gravitational experiments 
are still described by ordinary matter dominated theory and the PPN constraints must 
be satisfied. 



3 Solutions 



The general solution of Eq. ([15]) reads 



±x(t) 



-Cit 



(23) 



where M\ and M 2 are arbitrary constants of integration. These solutions fall into three 
classes, depending on the sign of constant C = C\ + 2C 2 . A detailed classification and 
characterization of the solutions, including the phase portraits was given in Ref . [1] , here 
we focus only on these solutions which approach the GR limit asymptotically in time. 



3.1 "Exponential solutions" 

In the case C > solutions read 



± x = e 



-dt 
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If C\ > and C 2 < the solutions exponentially converge to the GR limit, behaving as 



± x\ 
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Here we have denoted the constant of integration as M 1 = e 2* 1 (v / <5 Cl ) for some arbitrary 
moment t\. All solutions satisfy an asymptotic condition 
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The Hubble parameter reads 
H 



q 1 e -(Ci-V5)(t-ti) 
^ = — ± — 
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and effective barotropic index ([15]) can be calculated 



•1 ± 
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Now we can determine whether a model in the theory characterized by distinct param- 
eters (Ci,C = C\ + 26*2,^4^) approaches the de Sitter spacetime from quintessence side 
(w e ff > —1) or from the phantom side (w e g < —1). Approximate expressions of the 
PPN parameters f f20|) . ( 12 ip indicate that they approach the GR values /3 = 1, 7 = 1 
exponentially 

Solutions ( I24p may have interesting features at certain finite moments of time. Firstly, 
if the theory allows phantom divide, i.e. if K p m in Eq. (Tl9|) are real numbers, then 
solutions ( f2"4j) may cross the phantom divide no more than at two moments t±: 



In 



M 2 (C 1 + K ± + VC) 
M 1 (C 1 + K ± - 



C) 



(29) 



Secondly, at finite moments U some solutions can achieve x(tb) = 0, x{tb) = depending 
on values of integration constants Mi, M 2 . The latter ones can be given in terms of initial 
values e.g. at t — 0, x(0) = xo, x(0) = x$. Indeed, there is a moment U such that 
x(tb) = 0, x(tb) = 0, namely 



1 , M 2 1 

t6 = c ln M7 = c ln 



X + Xq(Ci - C) 

xo + xoid + C)^ 



(30) 



Phase trajectories have a vertical slope there and can be described as "turning back" if 
we consider solutions ( |24|) with only one sign (+ or -). The region of initial values for 
trajectories with this property can be found from the condition 



x + x (Ci - C) 
Xo + Xoid + C) 



> 0. 



(31) 



On the half-plane x > it is a region between two rays Xq = Xq(C — C\) and Xq = 
—Xq(C + Ci). The corresponding region on the other half-plane x < is its mirror image 
with respect to the origin x = 0, x = 0. 



3.2 "Linear exponential solutions" 

In the case C = the solutions read 



± x 



-dt 



e^ Clt H - M 2 
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with Mi = ez Cltl . If C\ > the solutions exponentially converge to the GR limit, 
behaving as 

±x\ t ^ OQ =t 2 e- c ^\ (33) 
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the Hubble parameter reads 



1 3 ^ + 



and w e fj behaves as 
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However, this case is rather finetuned by the condition C\ = — 2C 2 . 



(34) 
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3.3 "Oscillating solutions" 

In the case C < the solutions read 



± x 



-Cit 



N, sin(~t^C[) - N 2 CO s(h^/\C\) 



(36) 



where Ni, N 2 are integration constants. In terms of the phase space (x,x) solutions do 
not have a definite slope at approaching asymptotically (t — > oo) to x = 0, x = 0. As 
they spiral to it at t — > oo their phase trajectories cross this point infinitely many times 
at moments t n , 



tr, 



N 2 

arctan — + nrc 
Ni 
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The spiral, however, must lie in one half-plane of domain of definition, either x > or 
x < 0. Approximate expressions of the PPN parameters ( 1201) . ( 12~TT) now reveal damped 
oscillatory behaviour around the GR values. 

If the theory allows crossing the phantom divide, i.e. if K± in Eq. ([191) are real, then 
the possible moments t± n of crossing occur on each winding of the spiral: 



\c\ 



+ (C 1 + K ± )N 2 + ^\C\N l 
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If K± is imaginary, then the effective barotropic index w e g stays below or above — 1. 



3.4 A hint for the full theory 

We have considered solutions ( |23l) of an approximate equation (1151) for the scalar field. 
Among them the most peculiar seem to be "oscillating solutions" and a question arises 
whether this type of solutions really occur in the full theory (jSJ)-©, or perhaps they 
are only an artefact of approximation. In our previous paper we argued on the basis 
of the phase space that the topology of trajectories in the nonlinear approximation is 
representative of those of the full system, and thus the correspondence holds pQ. 

Here we supplement this discussion by a specific yet rather unphysical model, where 
an additional indication can be found that such solutions indeed exist in the full theory. 
Let us choose the two functional degrees of freedom to be 

3\j> 

= 2(1 - fr) ' kV(*) = o(1-*), a = const>0. (39) 
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In the corresponding approximate theory we have \I/ G (0, 1], = 1, A* = — ~, Ci = 0, 
C 2 = — y> C = — T' an d thus the approximate solutions should belong to the oscillatory 
type. 

In the exact theory the Friedmann constraint ([3]) reads 



2 rr* * a(l-tf) 



H 2 + H ; - -i '- = . (40) 

Considering a subspace H = const in the phase space (H, \I>, Eq. ( 1401) is a second 
order algebraic equation for the phase trajectories 

(#2 + ^ ^2 + + 1^2 _ ^2 + ^ ^ _ ^ + « = q (41) 

and their topology (ellipses or hyperbolas) is given by the coefficients of the second order 
terms: 

H 2 -A(H 2 + -)- = --<0. (42) 
v 3 ; 4 3 V ' 

It follows that the phase trajectories must indeed be ellipses, i.e. we have oscillations. 



4 An example 

To illustrate the behaviour of oscillating solutions, let us specify the coupling function 
and potential as 

= ^r^y > kV W = ■ ( 43 ) 

Here ^ G (0, 1] and the 'GR point' is at = 1. It makes sense to measure the oscillations 
in the units of the analogue of Hubble time, T = H+t = Q-t. Using prime to denote 
derivation with respect to T, the integration constants iV\ and N2 are related to the 
initial conditions of xq = x(0) and x' = x'(0) as 

m= ^ S- , AT 2 = v /3^. (44) 

To satisfy the current weak field constraints f|T0|) - f|T2|) we choose the initial conditions 

Xq = -0.000025 , x' = 0.0025 . (45) 

Substituting the solution ( 1361) and constants ( H4l) . ( 1451) into the expansion ( lT8i) we obtain 
the evolution of w e ff, which is plotted on Fig. 1. Here we see the effective barotropic index 
going through the phantom divide line in a damped oscillatory manner. We also note that 
the period of oscillations is comparable to the Hubble time, roughly the same order of 
magnitude as the age of the Universe. This suggests a possible scenario for the universe, 
where during the matter domination era the scalar field has already relaxed close to the 
GR limit [12] , while later when the cosmological energy density of the potential becomes 
more significant, we should observe the effects of slowly oscillating dark energy. 
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Figure 1: An example of oscillating dark energy described in Sec. 4. 

5 Concluding remarks 

We have proposed and solved a nonlinear approximate equation (|15|) for decoupled scalar 
field in the framework of STG FLRW cosmological models describing the era when the 
energy density of the scalar potential dominates over the energy density of the ordinary 
matter and the universe has evolved close to the limit of general relativity. The behaviour 
of solutions which approach general relativity can be classified under two characteristic 
types: (i) "exponential" or "linear exponential" convergence, and (ii) damped "oscilla- 
tions" around general relativity. These solutions also determine the evolution of cosmo- 
logical and PPN parameters. Depending on the model, exponential solutions may cross 
the phantom divide line at most twice before approaching w e g- = — 1 from either above or 
below. In the oscillating type the dark energy effective barotropic index oscillates either 
in the quintessence regime, phantom regime, or crossing the phantom divide line once or 
twice during each period. The latter case is illustrated by an explicit example. 
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